///&/! 


NASA Contractor Report 178010 nasa-cr-178010 

19860006158 

ICASE REPORT NO. 85-47 


ICASE 


TAYLOR— GORTLER INSTABILITIES OF TOLLMI EN- S CHLICHT I NG WAVES 

AND OTHER FLOWS GOVERNED BY 

THE INTERACTIVE BOUNDARY LAYER EQUATIONS 


Philip Hall 


Y. i 1 i 

u tii 



t 5 


Contract No* NAS1-17070 JAN 6 198 j 

October 1985 

LANGLEY RESEARCH CENTER 
LIBRARY, NASA 
HAMPTON,. VIRGINIA 


INSTITUTE FOR COMPUTER APPLICATIONS IN SCIENCE AND ENGINEERING 
NASA Langley Research Center, Hampton, Virginia 23665 

Operated by the Universities Space Research Association 


(W\SA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton. Virginia 23665 



3 1176 00189 9641 



TAYLOR-GORTLER INSTABILITIES OF TOLLMIEN-SCHLICHTING WAVES 


AND OTHER FLOWS GOVERNED BY 
THE INTERACTIVE BOUNDARY LAYER EQUATIONS 


Philip Hall 

University of Exeter, Exeter, England 
and 

Institute for Computer Applications in Science and Engineering 


ABSTRACT 

The Taylor-Gortler vortex instability equations are formulated for steady 
and unsteady interacting boundary layer flows of the type which arise in 
triple-deck theory. The effective Gortler number is shown to be a function of 
the wall shape in the boundary layer and the possibility of both steady and 
unsteady Taylor-Gortler modes exists. As an example the steady flow in a 
symmetrically constricted channel is considered and it is shown that unstable 
Gortler vortices exist before the boundary layers at the wall develop the 
Goldstein singularity discussed by Smith and Daniels (1981). As an example of 
an unsteady spatially varying basic state we also consider the instability of 
high frequency large amplitude Tollmien-Schlichting waves in a curved 
channel. It is shown that they are unstable in the first "Stokes layer stage" 
of the hierarchy of nonlinear states discussed by Smith and Burggraf (1985). 
The Tollmien-Schlichting waves are shown to be unstable in the presence of 
both convex and concave curvature. 
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INTRODUCTION 


Our concern is with the Taylor-Gortler instability of interactive boundary 
layer flows of the type which occur in triple-deck theory. Thus, we 

investigate the instability of the ^lower deck^ boundary layer which is set up 
when a classical boundary layer with Reynolds number Rg encounters a hump of 

c /o —3/8 

height and length order R^ and R^ respectively. We find that the 

form of the Taylor-Gortler instability equations in the lower deck are almost 
identical to those appropriate to a classical boundary layer. The main 
difference is that the wall shape function f(X,T) enters the instability 
equations and, in fact, for steady flows f^ plays the role of the Gortler 
number. 

We shall see that for unsteady interactive boundary layer flows both 
steady and unsteady Taylor-Gortler vortices of the type discussed by Hall 
(1982, 1983) and Serainara and Hall (1976) respectively are possible. These 
flows are also potentially unstable to short wavelength Rayleigh modes and the 
reader is referred to the papers by Bodonyi and Smith (1985) and Tutty and 
Cowley (1985) for a discussion of that problem. In general the instability 
equations which we derive must be solved numerically because a parallel flow 
and/or a quasi-steady approximation cannot be justified. However, in order to 
demonstrate that some of these flows are unstable we shall here concentrate on 
two problems for which some asymptotic progress is possible. 

Firstly, we look in detail at the steady flow in a symmetrically 

constricted channel. Smith and Daniels (1981) have shown that when h the 

scaled height of constriction becomes large a classical boundary layer of 
-l / 9 

thickness 0(h * ) is set up within the wall boundary layer. This inner 

boundary layer develops a Goldstein singularity beyond the minimum channel 
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width position, but Smith and Daniels showed that the singularity could be 
removed without any upstream influence being set up. Here we consider the 
instability of the flow before the boundary layer develops the singularity. 
For large values of h we are able to solve the instability equations 
asymptotically and demonstrate the instability of the h“ ^2 layer. In 
general the instability occurs before the Goldstein singularity develops; 
however, it is possible, as in Smith and. Daniels (1981), to choose humps with 
the required concave curvature only beyond the position where the singularity 
develops. 

Secondly, we look at the unsteady interactive boundary layer which governs 
the growth of Tollraien-Schlichting waves in parallel or boundary layer 
flows. Here the unsteadiness is characterized by £2 the frequency of the 
Tollmien-Schlichting wave. For definiteness, and to avoid the complications 
of boundary layer growth, we look at the instability of the waves in a 
slightly curved channel. Recently Smith and Burggraf (1985) looked at the 
structure of high frequency large amplitude Tollmien-Schlichting waves in a 
variety of situations. Dependent on the size of the disturbance and the 
particular flow under investigation, they found a hierarchy of nonlinear 
partial differential systems to describe the disturbance. The first nonlinear 
stage discussed by Smith and Burggraf is such that the disturbance has the 
form of a Stokes layer near the wall. We investigate the instability of this 
flow and identify the critical disturbance size above which the Tollmien- 
Schlichting wave is unstable to the Stokes layer Taylor-Gortler mode 
identified by Seminara and Hall (1976). When the Tollmien-Schlichting wave 
has amplitude greater than this critical value, a three-dimensional flow 
containing streamwise vortices develops; the effect of this new flow on the 
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growth of the wave into the larger amplitude states of Smith and Burggraf is 
beyond the scope of this paper. The procedure adopted in the rest of this 

paper is as follows: in Section 2 we derive the equations governing the 

« 

centrifugal instability of 'lower-deck' boundary layer flows. In Section 3 we 
show how these equations can be solved for the steady symmetric channel flow 
problem. In Section 4 we carry out a similar analysis for Tollmien- 
Schlichting waves in curved channel flows whilst in Section 5 we draw some 
conclusions. 


2. THE TAYLOR-GORTLER INSTABILITY EQUATIONS FOR TRIPLE-DECK FLOWS 

It is useful at this stage to discuss briefly the Taylor-Gortler 
instability equations for a classical two-dimensional boundary layer flow over 
a curved wall. The reader is referred to the papers by, for example, Gortler 
(1940), Smith (1955), Floryan and Saric (1979), Hall (1982) for a discussion 
of the approximations required to obtain a self-consistent set of linear 
stability equations. 

Suppose then that i and Uq are typical length and velocity scales for 
the flow and that v is the kinematic viscosity of the fluid. We define a 
Reynolds number Rg by 


“e = 


Un l 


( 2 . 1 ) 


and take x and y to be dimensionless variables measuring distance along 

and normal to a surface with local curvature a ^ K(x). The variable x is 

- 1 /, 

scaled on i whilst y is a boundary layer variable scaled on £R^ z . The 
basic flow u 0 is of the form 
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_ 1 u — 

= U 0 l u(x>y)>R E v(x,y),0) + 
and this flow is perturbed by writing 

~ = ~B + U o( U(x »y ) » R E 1/2v(x »y).R^ 1/2 W(x,y))exp{iR E 1/2 kz}. (2.2) 

Here k is a nondiraensional wavenumber in the spanwise direction, and we have 
assumed that the instability occurs on the boundary layer thickness length 
scale. From the momentum and continuity equations we can show that in the 

limit R E •> 00 with the Gortler number G = held fixed the linear 

stability equations are 

U x + V y + ikW = °» (2.3a) 

uU x + Uu x + vU y + Vu y = {3 y - k 2 }U, (2.3b) 

uV x + U \ + ^ V y + V ^y + K(x)G uU = -p y + {9y ~ k 2 }V, (2.3c) 

uW^ + vW y = -ikP + {3 2 - k 2 }W. (2.3d) 

Here P is the nondimensional pressure perturbation corresponding to (U, V. 
W) and we have assumed that the perturbation is steady. The generalization 
of (2.3) to a weakly three-dimensional boundary layer is given by Hall (1985). 

The essential difficulty with (2.3) is that for G and k 0(1) there is 
no rational reason why a parallel flow approximation should be made and the 
P a rtial differential system must be solved numerically as was done by Hall 
(1983). For k » 1, but G ~ an asymptotic solution to (2.3) was given 



by Hall (1982) who showed that in this limit nonparallel effects can be taken 
care of in a self-consistent manner. For 0(1) wavenumbers the numerical 
calculations of Hall (1983) showed that the position of neutral stability is a 
function of the initial disturbance. At higher Gortler numbers the local 
growth rate approaches the asymptotic result which, in this regime, is 
consistent with a parallel flow theory calculation. Unf ortunateley , it has 
been assumed elsewhere that this latter result justifies the use of parallel 
flow theories. However, in the only regime where the parallel flow theories 
are valid, i.e., k » 1, an asymptotic result of at least the same accuracy as 
any parallel flow theory can be written down in closed form with little 
effort. 

We now show how the equations corresponding to (2.3) can be derived for a 
basic flow governed by some interactive boundary layer structure. For 
definiteness we focus on a flow governed by triple-deck theory; the 
formulation for other structures is essentially identical. Consider then the 
flow over the wall 

y = e 5 f(X,T), 


where X = e x, and 
—1 /8 

parameter e = R_, 

Ej 


T = e 


-2 tU 0 


Here t denotes time whilst the small 


We define the lower deck variable Y by 


Y 


-5 

e 


y> 


and in the lower deck the basic state expands as 


-^BL 


= (eu(X,Y),e 3 v(X,Y),0) + 

0 
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whilst the pressure expands as 


pi. 2 — . . 

^j-= e p(X) + .... 


The equations which determine the flow in the lower deck 


are 


U T + u U X + V u Y = - p x + u YY , 


(2.4a) 


u x + V Y = °> 


(2.4b) 


whilst the boundary conditions at the wall are 


u = 0, v = f on Y = f(X,T), 


(2.5a) 


and at infinity we require 


u ->■ Y + A(X,T), 


(2.5b) 


where A is the displacement function which must be related to p through a 
pressure - displacement law. If we make the unsteady Prandtl transform 


Y Y + f(X,t), v -*■ v + uf + f , 

A J. 


( 2 . 6 ) 


then (2.4a), (2.4b) are unchanged whilst (2.5) reduces to 


u = v = 0, Y = 0, 
u-vY+A + f, Y + co. 


(2.7) 
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We now investigate the instabilty of this flow to Taylor-Gortler vortices 
which might be associated with either the steady or unsteady components of the 
basic state. We look for perturbations confined to the lower deck having 
spanwise wavelength comparable with the lower deck thickness. The possible 
source of instability is, of course, the curvature of the wall in the lower 
deck. We write 


^-= u BL + A(eU(X,Y,T), e 3 V(X,Y,T), e 3 W(X,Y,T))E, 


( 2 . 8 ) 


where A <X 1 and E = exp — — . Here Z has been scaled on £Rg ^ and we 

e 

have assumed in (2.8) that the normal and spanwise velocity components are 
comparable. This is the usual case for Taylor-Gortler instabilities, and the 
relative scaling of the X and Y velocity components is again consistent 
with that appropriate to the small gap Taylor vortex problem (see, for 
example, Davey (1962)). The pressure perturbation p + in the lower deck 
expands as 


+ 

p a 



e 6 P(X,Y,T)E, 


(2.9) 


and the above relatively small scaling for P enables us to retain the 
convective and diffusive terms in the Y and Z momentum equations. It 

remains for us to substitute the perturbed flow into the Navier-Stokes 
equations for the lower deck and with A « 1 linearize about the basic 
state. We note that at this stage it has not been necessary to define a 
Gortler number for the flow. After linearizing about the basic state and 
making the Prandtl transform (2.6) together with 
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V ^ V + uf x , 

we obtain 

U x + V y + ikW = 0, 

U T + uU X + Uu X + vU Y + % = {9y “ k2 } U » 

V T + uV x + Uv x + vV Y + v7 y + 2U{f xx u + f XT } = -P Y + f 9 2 - k 2 } V , 

W T + uW X + ^ W Y = “ ikP + 1 9 Y ~ k2 J W> (2.10) 

which must be solved subject to 

u = v = W = 0, Y = 0, < 0 . (2.11) 

We see that the generalization of (2.3) to an unsteady triple-deck flow leads 
to almost the same equations but with KGuU replaced by 2{f u + f }. For 
steady triple-deck flows this means that f xx plays the role of the Gortler 

number, whilst for unsteady flows an extra term proportional to Uf^ arises 

due to the vertical motion of the boundary. For time-periodic basic states 
the system (2,10) contains the terms which lead to centrifugal instabilities 
in Stokes layers so in general we must be alert to the possibility of both 
types of Gortler instability. 

The solution of (2.10) is in general a numerical problem which, for time- 
dependent flows, will be an order of magnitude more difficult than the steady 
state calculations of Hall (1983). In the next section we will look at a 

particular steady state which it is possible for us to solve (2.10) 



-9- 


asymptotically in a self-consistent manner. Finally we note that, with the 
appropriate scalings used, (2.10) apply to other interactive boundary layer 
flows governed by (2.4), (2.7). 

3* SYMMETRIC CHANNEL FLOWS 

In general it is not possible to solve (2.10) analytically; however, we 
now show how asymptotic methods can be used for a particular steady basic 
flow. We refer to the internal channel flows discussed by Smith and Daniels 
(1981). In that problem the wall boundary layer thickness is ~ R ' and 
the x variations are on an 0(1) length scale. However, if the disturbance 
quantities are scaled appropriately then (2.10) still apply. The basic state 
satisfies (2.4) but with A = 0 in (2.7) since the channel is symmetric. 

For steady flows we saw in Section 2 that f xx plays the role of the 
Gortler number so that we expect the flow to become more unstable with 
increasing hump height. In view of the work of Hall (1982) we might then 
expect that an asymptotic solution of (2.10) should be possible. The Smith- 
Daniels problem is a suitable candidate for such an analysis because its 
structure for |f| » 1 is reasonably well understood. Suppose that we write 

f(x) = hF(X), (3.1) 

with h » 1. In this situation a classical boundary layer extending to 
X = -to and of thickness ~ h ^ is attached to the hump. In this layer 
u, v, expand as 


u = hU + 


(3.1a) 
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v 



(3.1b) 


where U and V are functions of X and £ = h Y , and satisfy the 

classical boundary layer equations with pressure gradient -FF X and TT + F 

when £-»■<». The effective Gortler number for the inner 0(h~^ 2 ) boundary 

7 /2 

layer then becomes 0(h ) so that, on the basis of the work of Hall (1983), 

we expect that neutral modes will have k ~ 0(h^ 8 ). Moreover, these modes 

will be concentrated in an internal layer of depth which is located 

so as to maximize the local spatial growth rate. At the location X we 

assume the layer is centered on £ = Y( x ) and write 


q = h 3/,16 (? - £(X)). 


(3.2) 


The wavenumber k is then expanded as 

k = kg Xi ^ + kj h ^ + •••, 


(3.3) 


whilst we write 


U = {U 0 (q,X) + h 3/16 UjOl.X) + h 3/8 u 2 (n,x) + • • • 


(3.4) 


together with similar expansion for 
& 

quantity E is defined by 


V W 

h U4 • „1/16 


and 


h i5/i6 • 


Here the 


* 

E 


exp 



X 

/ {e 0 (x) + h 


-3/8 


0 j(X) + ...}dx 




(3.5) 
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so that { B (X) } determine the spatial growth of the disturbance. In fact, 
we will concentrate on the neutral case and to the order which we proceed here 
it is not necessary to distinguish between the growth rates for different flow 
quantities. Finally, near 5 the basic state expands as 


u = h{U Q (X) + h 3/16 nU 1 (X) + h -3/8 n 2 U 2 (X) + •••} (3.6a) 

V = h 1 / 2 {V Q (X) + h -3/16 nV^X) + h -3/8 n 2 V 2 (X) + •••} (3.6b) 


where 


U.(X) = 

1 1 ! 


v t (x) - v -l p s l 


(3.6c) 


It remains for us to substitute the above expansion into (2.10) and 

“3/16 

successively equate like powers of order h • The zeroth order problem is 


found to be 


u o + k o) v o + 2F xx u o V °« 


(3.7a) 


0*0 U 0 + k 0) U 0 + V 0 U 1 = °’ 


(3.7b) 


iP 0 + k 0 W 0 " °> 


(3.7c) 


v 0 + ik 0 w 0 = °* 


(3.7d) 


The required consistency condition for (3.7a), (3.7b) yields the zeroth order 
eigenrelation 

(®o u o + k o^ = u i u o F xx» 


(3.8) 
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and the potentially unstable root of this equation is 


u o 6 o " k o “ ^ u i u o F xx » 


which exists only for U, U n F YY > 0. At a position where U, U, F 

i u aa 1 0 XX 

vanishes we have a coalescence of modes and a transition region is required. 
The most unstable position in the layer is such that 


28 0 U 0 (6 0 U 0 + k£) = [Uj + 2U„ 0 2 ] Fxx , 


(3.9) 


and, with f, k Q given, (3.8) , (3.9) fix c and 0 . 

restrict our attention to the neutral case and set 8^ = 6 

0 P 1 

the flow is neutral at X if 


k 0 ■ U 1 U 0 F XX’ 


From now on we 
$2 = 0 so that 


(3.10) 


and £ is determined by the condition 


u ? + 


2U 0 U 2 - 0, 


which requires that |u IK | has a maximum at ? = At this stage the 
vertical structure of the eigenfunction is not determined; at second-order the 
following equation for Vq emerges as a solvability condition: 


~ 3 k 0 k l V 0 + k (T ^ F Xxt U l U 2 + U 3 U 0 ] V 0 = °‘ 


-2 2 


(3.11) 


The solutions of (3.11) which decay when ri 


-»■ oo are 
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v 0 = V Qn = U(-n -V 2 > YH> (3.12) 

where U(-n - V 2 » yn) is a parabolic cylinder function and 

„ - /21-FjUj 0, + U 3 FJk- 2 } ^ , 

the wavenumber k, must then satisfy 

k l = k ln = Ikjj’ ( n +1/ 2 1* (3.13) 

1 /.-> 

The most unstable mode corresponds to n - 0 so that correct to order h /z 
the neutral wavenumber is 

gi "2 * “pi 

"0 “1 

Thus for a given hump size h and shape f the flow is neutrally stable at 

X to a vortex flow with wavenumber k if (3.14) is satisfied. Clearly 

instability occurs only for F^. > 0, and if we restrict our attention 

to humps with f > 0 then upstream separation does not occur and it can be 
shown that until the Goldstein singularity develops IJ IJ > 0 so that only a 
locally concave wall can lead to instability. 

If we assume that F(X) and k » 1 are given, then an alternative 
interpretation of (3.14) is necessary. In this situation we can think of 

(3.14) as an implicit equation for the hump height h » 1 which makes the 

Gortler vortex flow with wavenumber k » 1 neutral at some position in the 
flow field. In the following discussion we assume that F 0, X + and 


k ■ ( U 1 u o F xx) 1/4 1,7/8 + T 


I >1 


V2, 

T • • • . 


(3.14) 
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that F xx is positive in (-», -c) and negative in (-», d) where C and 
D are positive constants. Since we also know that 

u = 0, ? = 0, 0, c + - 

and 

U ~ 0(f), X ^ 

it follows that F xx has at least one maximum in < x < -C, 

0 < £ ^ "• Suppose that the largest maximum occurs at 


x ■ V ' - V 

then ignoring the order h ^ term in (3# 14) we see that the minimum hump 

height h^ which leads, to a neutral vortex anywhere (in fact at (X r , r )) 
is given by 

h C " <X C (lt C )F Xx )" 2/7 k 8/7 . (3.15) 

If h is increased beyond h^ there will be two neutral locations at 
~ C ~ a > -c + S with a, 0 > 0 each corresponding to the wavenumber k. 
Between these positions the flow is formally unstable; the instability will 
amplify by an amount 0(exp[h 3/4 I]) for some I > 0 in this interval and 
will become nonlinear if its initial size is sufficiently large* For humps 
with h ~ 0(1) we expect that a similar situation arises for k - 0(1) with 
at least one finite interval where instability occurs* Beyond the position of 
the maximum constriction a further region where instability can occur will 
exist so long as the Goldstein singularity is not encountered* 



4. THE INSTABILITY OF TOLLMIEN-SCHLICHTING WAVES IN CURVED CHANNEL FLOWS 


In the previous section we described how a particular steady solution of 
(2.4), (2.7) becomes unstable to steady Taylor-Gortler vortices. We now show 
how a time-periodic solution of that system can also become unstable to a 
Stokes layer Taylor-Gortler vortex. The particular type of time-periodic 
basic state which we consider corresponds to a large amplitude high frequency 
Tollraien-Schlichting wave propagating in a curved channel. Here the curvature 
which causes the instability is not on the triple deck length scale, so it is 
necessary to say a few words about the derivation of the appropriate form of 
the equations corresponding to (2.10). 

We consider the flow driven by an azimuthal pressure gradient between 
concentric cylinders of radii a, a + d. The maximum flow velocity is taken 


to be 



and if we define dimensionless variables 


x and y by 


x 


a0 

d 



0 , 


y 



- d 


then in the absence of either Tollmien-Schlichting waves or Taylor-Gortler 
vortices the basic state for 6 « 1 is 

Ug = U 0 (u Q + 0(6), 0, 0) (4.1) 

with 

Uq = y(i - y). (4.2) 


The Reynolds numer Rg is defined by 
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and Dean (1928) showed that (4.1) is centrifugally unstable for 6 ~ 0(1). 
The flow is also unstable to Tollraien-Schlichting waves for ^ > 5774. We 

ignore the steady Taylor-Gortler mode and examine the instability of finite 
amplitude Tollmien-Schlichting waves to Stokes layer Taylor-Gortler vortex 
modes. These occur for 6 ~ 0(1) and so are apparently less important 

than the steady vortex mode. However, we shall see that they occur both near 
the inner and outer cylinders so that in external flows over convex walls we 
can expect this mode to be the only centrifugal one available. Our choice of 
the curved channel flow problem rather than an external boundary layer enables 
us to investigate the instability mechanism without the possibly insurmount- 
able difficulties of handling the effect of boundary layer growth. 

We first describe the large amplitude high frequency Tollmien-Schlichting 
disturbances to (4.1) which exist for Rg » 1. The discussion we give is 
taken from the recent work of Smith and Burggraf (1985) who investigated a 
range of flow regimes where such disturbances can exist. We take z and t 
to be dimensionless axial and time variables scaled on d and U 0 / d . The 
appropriate length scale for a Tollraien-Schlichting wave in a channel is 
0(R^7) so that if 

we define 



and 

X = ex (4.4) 

then the wall layers at y = 0,1 are of thickness 0(e 2 ) so that near 
y = 0 we write 


Y = y/e 2 . 


(4.5) 
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Following Smith and Burggraf we seek a flow for Y ~ 0(1) of the form 


u = U Q (e 2 u(X,Y,T), e 5 v(X,Y,T),0) + 


(4.6) 


where T = e t. The corresponding pressure perturbation is 
and u, v, p are determined by 



4 

e 


p(X,T) 


u T + u u x + v u Y = -P x + u YY 

^X + v Y = ° 

u = v = 0, Y = 0 


u Y + A(X,T) , Y oo 


(4.7) 


Near y = 1 a similar boundary layer exists whilst in the core we have 


U = U o( u o ,0, °) + U 0 ( e2 6)3 v(X,y,T),0) + 


P 


” U 0 4 

d £ 


p(x,y,T) 


+ • • • # 


(4.8) 


Here u, v, p satisfy 


u o u x 





°X + v y 


= 0 


u o v x = -Po 


(4.9) 
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and the appropriate solution is 


(u,v) = (A(X,T)u 0y , - A^X.Duq), 


(4.10) 


Po 



(4.11) 


so that the required pressure-displacement relationship is 


y=l 


" Pr 


1 y=-0 


30 A XX* 


(4.12) 


We are interested in the large amplitude high frequency solutions of (4.7) 
discussed by Smith and Burggraf (1985). The latter authors investigated a 
hierarchy of high frequency large amplitude states beginning with the case 
■gY ~ 0(52 ) » 1, Yx ~ 0(ft 2), u ~ 0(1). We shall concern ourselves here only 
with the latter state and the reader is referred to the Smith-Burggraf paper 
for a discussion of the remarkable range of more nonlinear states which occur 
for u » 1. 

For ft » 1 we write 


n 


= ft 


l h 


Y, 


3^, ft3 + ft ^ 3 + ••• 

T T 

V? 

3 V -»■ ft 11 3 

X X 

and note that linear stability would then lead to instability on the T 
scale. We note that Smith and Burggraf looked for spatial rather than 
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temporal growth. Here it is convenient to look for growth in time because it 
leads to an eigenvlaue problem already solved in a different context. The 
spatially growing case, except in the locally neutral state, leads to a 
different eigenvalue problem. The work of Smith and Burggraf shows directly 
that in this regime 


where 


u = Uq^ + 0(ft ^ 2 ) 

(4.13a) 

v = v qi + o(cf 1/2 ) 

(4.13b) 

= {p 01 (l - e mn )E + C.C.}, 

(4.13c) 

p mn 1 

- (icpoiCi - S + „) E + C - C -}- 

(4.13d) 


Here C.C. deontes "complex conjugate" and E denotes expi[a X - T} with 
a given by the eigenrelation 

a = /30, 

3 fT i 

whilst m = e — . 

4 

Thus the Tollmien-Schlichting wave is bigger than the steady basic state 
for r\ ~ 0(1) and is just a Stokes layer flow. At higher order the Smith- 
Burggraf approach shows that p Q1 grows exponentially with T and does not 
equilibrate. We will look at the instability of the Tollmien-Schlichting wave 
during this stage of its growth. 

We shall firstly consider the instability of (4.6) to a Taylor-Gortler 

2 

vortex perturbation with axial wavelength 0(e ). We thus write 

u = U 0 (e 2 u + e 2 [U(X,Y,T)exp ikze” 2 + C.C.] + •••), (4.14a) 
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v = u 0 ( e5 v + e 5 [V(X,Y,T)exp ikze“ 2 + C.C.] + ...), (4.14b) 

5 —2 

w = U o( £ [W(X,Y,T)exp ikze + C.C.] + •••), (4.14c) 

P = (e 4 P + e 10 P(X,Y,T)exp ikze“ 2 + C.C.] + (4.14d) 


and then 

A 

26 = De 

where D can be interpreted as a Taylor-Gortler number for the Tollmien- 
Schlichting wave. The equations (4.14a), (4.14b), (4.14c), (4.14d) are then 
substituted into the momentum and continuity equations, and after linearizing 
about the basic state we find that (U, V, W, P) satisfies (2.10) but with 

the term 2U{f^ u + f^} replaced by Dull. At this stage the eigenvalue 

problem D — D(k) could in principle be solved numerically for any given 

basic state. We shall proceed by looking at the high frequency limit of the 
Tollmien— Schlichting waves in order to make some analytical progress. We 

stress that we expect instability to occur for 0(1) frequencies but do not 
pursue the necessary large scale computational task required to verify this 
speculative remark. 

Let us now look at the high frequency limit of (2.10) and show that the 
eigenvalue problem governing the instability of the growing Tollmien- 
Schlichting waves can be reduced to one already studied in the context of 
Stokes layer instabilities. 
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Thus, in the Stokes layer we expand U, V, W, P in the form 

T 

U = [U Q (X,n,T) + 0(n - ^ 2 )]exp fi 3/ ^ 2 / cr(T)dT , (4.15a) 

V = [n 1/2 V 0 (X,n,T) + 0(n°)]exp fi 3/2 / o(T)dT, (4.15b) 

W = [fi 1/2 W 0 (X,n,T) + 0(n°) ]exp fi 3/2 / <j(T)dT, (4.15c) 

rsj 

T 

P = [S2P Q (X,n,T) + 0(S2^ 2 )]exp ft 3/2 / o(T)dT (4. 15d) 

whilst k and D expand as 

k = k Q Q 1/2 + 0(n°), (4. 15e) 

D = D q fi 3 ^ 2 + 0(a). (4 . 15f ) 

The X variation enters the zeroth order problem only through the X 
dependence of u 01* T ^ e S rowt ^ r ^te a(T) depends on T through the T 

dependence of Pq^. 

If the expansions (4.14), (4.15) are substituted into (2.10), the zeroth 
order problem is found to be 

v o„ + lk 0 “o = °> <4 - 16a) 

au 0 ♦ U _ + v 0 u 01n - (a 2 n - kg}U 0> 


(4.16b) 
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° v o + v oT - “o “01 D oi ■ - p o„ + I 3 n - k ot v o- (4a6c) 

CW 0 + H _ - -lt 0 P 0 + {9 2 n - k 2 }w 0 , (4. 16d) 

which must be solved subject to 


U 0 = V 0 = W 0 = °» (4.17a) 

Uq> ^q> w q + 0> n 00 (4.17b) 

so that the vortex structure is confined to the Stokes layer. The eigenvalue 
problem specified (4.16), (4.17) is identical to that studied by Hall (1984) 
in the context of Schlichting's (1932) transversely oscillating cylinder 
problem. In fact, we identify (k,T) of (2.12) of Hall (1984) with (/2k 0 , 

,/2D o I P 01 f )• 

It follows from the numerical calculations of Hall (1984) for the neutral 
case a=0 that the Tollmien-Schlichting wave is formallly unstable for 

D 0 I Pqi I 2 > 8.48 (4.18) 

which, for a given value of Dq, determines IPqiI the critical Tollmien- 
Schlichting wave amplitude. 

Not surprisingly a similar analysis governs the instability problem for 
the Tollmien-Schlichting wall layer at y = 1. The only change is that, 
because the layer is now on a concave wall, the sign of D in (4.16) must be 
switched. Papageorgiu (1985) has investigated that eigenvalue problem, and 
the critical state is then determined by 
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D p 01 T > 5.02. 


Thus the layer at y = 1 becomes unstable first, and presumably the flow 
becomes three-dimensional before the larger amplitude two-dimensonal states of 
Smith and Burggraf develop. It is known from the experimental work of 
Seminara and Hall (1976) and Park, Barenghi, and Donnelly (1980) that the 
initial Stokes layer instability is followed by a secondary mode of 
instability at about 30% above the first critical Taylor number. In this 
regime the vortices interact and the disturbance persists beyond the Stokes 
layer. No adequate theoretical description of this nonequilibrium state is 
yet available, but the consequences for the problem discussed here are 
important. We refer to the fact that if the Tollmien-Schlichting wave also 
undergoes this secondary mode of instability then there will be a mechanism 
for disturbances inside the Stokes layer to penetrate outside the boundary 
layer. 


We further note that if the mechanism described here does indeed apply to 
external flows then Tollmien-Schlichting waves will generate Taylor-Gortler 
vortices if convex or concave regions exist. In the concave regions the 
Taylor-Gortler mechanism associated with the main deck basic state will be 
more unstable so that the Tollmien-Schlichting breakdown into the Stokes layer 
mode is probably only of practical importance in convex regions. 

Finally, we note from the calculations of Hall (1984) that when 
instability occurs , Wq , Pq are of the form 




5. CONCLUSION 

We have shown that interactive boundary layer flows of the type which 
arise in triple-deck theory can support Taylor-Gortler vortices. The form of 
the equations for steady flows is identical to that found for classcial 
boundary layers if we interpret the wall curvature as the Gortler number. For 
unsteady boundary layers, an extra term proportional to the streamwise 
gradient of the wall velocity is introduced into the equations. 

We have seen that both steady and time— periodic boundary layer Gortler 
vortices can be described within the above framework. In particular, we 
showed that a large amplitude high frequency Tollmien— Schlichting wave can 
interact with a curved wall to give a Stokes layer Gortler vortex. If such a 
result also holds for external flows, it means that wall curvture, either 
concave or convex, can lead to the breakdown of two-dimensional Tollmien- 
Schlichting waves into a three-dimensional flow having streamwise vortices in 
a sublayer. 

In general, the partial differential equations describing the evolution of 
Taylor-Gortler vortices in interactive boundary layer flows will have to be 
solved numerically. The results we found for the Smith-Daniels flow in a 
symmetric channel suggest that for 0(1) values of h instability will 



-25- 


occur. Thus we expect that a numerical investigation of the instability 
equations for h ~ 0(1) will show the presence of growing Gortler vortices 
for some range of values of X whilst for large h we would expect to 
recover the asymptotic results of Section 3. 

It is not yet known what effect the above type of instability will have on 
separation; certainly linear stability theory can lead to no predictions on 
this matter. However, in situations where the small wavelength analysis of 
Hall (1982) can be used we know that the instability will be concentrated away 
from the wall. Thus for the Smith-Daniel problem we can show that when the 
Goldstein singularity develops the instability still persists since the new 
structure required to remove the singularity is confined to a region near the 
wall. 

We should not overlook the fact that the flows which we have investigated 
are also susceptible to Tollmien-Schlichting waves. In a situation where the 
basic flow can support both modes of instability we expect that, in view of 
the much larger spatial growth rates, the Gortler mode will dominate. An 
adequate description of even a weakly nonlinear interaction between Tollmien- 
Schlichting waves and Gortler vortices has not, to the author's knowledge, 
been given for any flow. 
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